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Abstract 

We have performed a systematic study of the phase transition in 
the pure compact U(l) lattice gauge theory in the extended coupling 
parameter space 7) on toroidal and spherical lattices. The obser- 
vation of a non-zero latent heat in both topologies for all investigated 
7 € [+0.2,-0.4], together with an exponent v e g ~ 1/d when large 
enough lattices are considered, lead us to conclude that the phase 
transition is first order. For negative 7, our results point to an in- 
creasingly weak first order transition as 7 is made more negative. 
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1 Introduction and overview 



In this work the 4D pure compact U(l) gauge theory is revisited. This 
model, being the prototype of gauge abelian interaction, has received a lot 
of attention through the years, since it constitutes a preliminary step in ap- 
proaching QED. However, the questions related with the mechanisms driving 
the transition and with the order of the phase transition itself, have turned 
out to be a more controversial issue than one might expect from its apparent 
simplicity. 

Through this section we shall summarize the main results obtained up to 
date. 

The 4D pure compact U(l) theory is known to undergo a phase transition 
between a region in which the photons are confined (strong coupling) and 
another region in which they carry a Coulomb potential (weak coupling). 
The theory can be discretized in the lattice using the Wilson action 

S w = (3j2l 1 -RetrC/p] . (1) 
p 

Most of the work has been done on the hypercubic lattice with periodic 
boundary conditions, namely on the hypertorus. On this topology the pioneer 
studies (see eg. |], ^|, [|) pointed to a continuous transition since numeri- 
cal simulations did not evidence signals of met ast ability. The confinement- 
deconfinement transition sets in at a coupling (5 ~ 1. The largest lattice size 
simulated in those works was L = 5. However, as simulations in larger lat- 
tices became available to computer resources, metastabilities were evidenced, 
and two-peaked energy histograms could be observed from L = 6 on (see eg. 
@j H) 11) pointing to a first order character for the phase transition. This 
idea was also supported by several works that approached the problem from 
the point of view of the Renormalization Group ]7|, j|, ||, [10| . 

Nowadays, though the general belief is that the transition is first order, 
there are still some puzzling questions. Namely, the numerical simulations 
up to L = 12 do not reveal a stable latent heat, it is rather decreasing 
with the lattice size, and also, the associated critical exponent v is always 
in the interval v = 0.29 — 0.32, which is different from the first order value 
[y = 0.25) and from the trivial second order one, namely v = 0.5. These 
facts led some authors to wonder about the possibility of a slowly vanishing 
latent heat, which could eventually go to zero in the thermodynamic limit. 
In short, the 4D pure compact U(l) theory could be critical provided the 
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two peaks appearing in the simulations are a finite size effect. Due to the 
physical implications of a critical theory with non trivial critical exponents 
this possibility was not completely ruled out. 

Now, we go back in time to introduce the extended Wilson action. Besides 
the Wilson action, an action including a term proportional to the square of 
the plaquettes was proposed by G. Bhanot at the beginning of the eighties 


Sew = P Y)- 1 ~ Re tr U ^ + 7 ~ Re tr ^pl ' (2) 
p p 

In that work the extended (/3/y) parameter space was explored in L = 4, 5, 
and it was suggested that the phase transition line /3 c (j) became first order 
for 7 > +0.2, while for 7 < 0.2 the transition appeared to be second order 
on those small lattices. The phase diagram is represented in Figure [I]. 

Simulations in larger hypercubic lattices have shown that two metastable 
states are present for all, up to date, positive and negative simulated 7 val- 
ues. What holds true is the observation that a positive 7 value reinforces the 
transition towards the first order, while a negative 7 value weakens the tran- 
sition. In fact, the latent heat, C\ at , decreases with 7, and then, a coupling 
where Ci at goes to zero was supossed to exist at some finite (negative) value 
of 7. This behavior induced some authors to conjecture about the existence 
of a tricritical point (TCP) at some negative coupling 7tcp where the order 
of the phase transition could change from being first order to be continuous 
P"2"fl . However, every simulation performed at negative 7 has exhibited dou- 
ble peak structures for large enough L. As for the Wilson action (7 = 0), for 
negative 7 a non stable latent heat is revealed by simulating in lattice sizes 
in the range L — 4, 12. 

From the point of view of the Renormalization Group this problem was 
studied in jT5] and the results pointed to a phase transition line /3 c (j) of first 
order for all 7 except at the end, 7 = — 00, where the transition is second 
order with trivial exponents. 

Nowadays, it is generally accepted that the transition is first order for 
positive 7 values. However for 7 < two scenarios arise: 1) The phase 
transition line (3 c {l) is first order everywhere except at 7 = —00; 2) There is 
a finite ^tcp < from which on, the phase transition is second order. 

As we have already remarked, all the previous discussions refer to simu- 
lations performed on the hypertorus. 

Recently some authors have cast a new hue on this problem. It was 
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suggested some years ago that the jump in the energy observed when the 
confining-deconfining transition takes place, could be due to the existence of 
closed monopole loops wrapping around the toroidal lattice, which due to the 



particular topology cannot be contracted to a point [14|. This hypothesis led 



authors to work on lattices homotopic to the sphere, since on this topology 



all monopole loops can be contracted to a single point [[15], [II], [LyJ . 

The relationship between the closed monopole loops and the disconti- 
nuity in the energy is far from being proved. The influence of those non 
trivial monopole loops on the phase transition has been studied by Rebbi 
and coworkers (18[. Their results suggest that the properties of the phase 
transition seem not to be affected by the existence of such closed monopole 
loops. 

However it is claimed the two state signal does dissapear when working 
on a lattice constructed considering the surface of a 5D cube, which has the 



topology of a sphere |T5|, |T7| . This single peak structure together with the 



estimation of the critical exponent v ~ 0.37 measured at 7 = 0, —0.2, —0.5, 
has induced authors in [H]] to claim that the double peak signal on the torus 
is a finite size effect, and when working on the sphere one observes that 
indeed the line /3 C (7) for 7 < is a line of second order phase transitions, all 
of them in the same universality class characterized by a v exponent around 
1/3. Assuming this second order scenario for the pure gauge theory, fermions 
have been included (on the toroidal topology) in order to study how fermions 
fit in the U(l) theory. Also in this second order phase transition is 

conjectured to exist |L9] . 



This situation is somehow disturbing, because one would expect to get 
the asymptotic behavior of a phase transition earlier on a lattice which is 
homogeneous and translationally invariant, namely the hypertorus, than on 
a lattice which is manifestly non homogeneous, as is the surface of a 5D 
cube, and on which the translational invariance is only recovered in the limit 
of infinite volume. 

First order phase transitions exhibit a wide variety of behaviors. Some 
of them do not show any pretransitional effect, the correlation length at the 
critical point, £ c , is small, and the thermodynamic quantities diverge as the 
volume for reasonable lattice sizes. However when the correlation length 
is very large, pretransitional effects take place, and they can also be very 
significant. In literature these transitions are referred to as Weak First Order 
(WFO) phase transitions and are characterized by small discontinuities, such 
SIS db small latent heat. 
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Depending on whether or not L is larger than £ c , a WFO phase transition 
presents two differentiate regions: 



• Transient region (L < £ c ). The transition behaves like a second order 
one, since £ c is effectively infinite. Thermodynamic quantities diverge 
less than L d and hence the effective critical exponent v is larger than 
1/d. 

• Asymptotic region (L ^> £ c ). The first order character of the transition 
is revealed by the onset of a stable two-peak structure and an effective 
v compatible with 1/d. 

How large should £ c be for the transition to be considered WFO seems 
a subjective matter. Since £ c 7^ always, all first order transitions would 
be weak in some extent. The key point for a transition to be considered 
WFO is the existence of pretransitional effects. The pretransitional behavior 
makes the transition along the transient region to look like a second order 
one with a definite critical coupling (3'. This second order phase transition 
is preempted by the true first order one at the coupling f3 c . The distance 
between the couplings f3' and (3 C can be taken as a measure of the weakness 
of the transition p0[ . Typical examples of such WFO transitions are 2D 
Potts models and the deconfining transition in finite T QCD ]2"T| . 

From this description it is clear that the observation of no double-peak 
structures does not prove by itself the second order character of a phase 
transition. 

In view of this, we decided to make a systematic study of the phase 
transition line /3 c (t) on the hyper-torus (HT topology hereafter) and on the 
lattice with spherical topology following |15], [IT] (which will be referred to 
as HS topology). On the hypertorus we improve on the statistics running 
up to L = 24 at some negative 7 values. On the sphere we have run lattices 
substantially larger than those employed by authors in [[HJ to check if the 
two-state signal reappears when larger spheres are considered. Our prelim- 
inary results, suggesting that the double peak structure indeed reappears, 



are published somewhere else [p3|| . Here, the statistics on those lattice sizes 
have been increased, and simulations in larger lattices are included. New 
observables are used to extract our conclusions. 
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2 Description of the model and observables 



We shall consider the parameter space described by the extended Wilson 
action, which can be expressed in terms of the plaquette angle in the following 

way: 

S = -/3^cosfl p -7^cos2# p . (3) 
p p 

We use for the simulations the conventional 4D hypercube with toroidal 
boundary conditions (hypertorus) , and, for comparison we also consider the 
surface of a 5D cube which is topologically equivalent to a 4D sphere. Con- 
trarily to the hypertorus, the surface of a 5D cube is not homogeneous. There 
are a number of sites which do not have the maximum connectivity, namely 
8 neighbors. Due to this fact, uncontrolled finite size effects are expected 
to turn up. Their influence can be somehow alleviated by the introduction 
of appropriate weight factors in those inhomogeneous sites [17]. Since the 
topology remains unchanged and we do not expect the order of the phase 
transition to be affected by the rounding, we do not use weight factors. 

As a notational remark, we label with L and N the side length for the 
HT topology and for the HS topology respectively. 

Next we define the energies associated to each term in the action 

iVp p 

E2 P = ^(J2 cos29 p) > (5) 



p p 



where N p stands for the number of plaquettes. On the hyper-torus this 
number is simply proportional to the volume, namely the forward plaquettes 
are N p = 6L 4 . On the sphere the number of plaquettes has a less simple 
expression, and can be computed as a function of N as N p = 60 (iV — l) 4 + 
20 (N — l) 2 . In this case, the system is not homogeneous and N p is not 
proportional to the number of points on the four dimensional surface, which 
is N 5 — (N — 2) 5 , some points having a number of surrounding plaquettes 
less than the possible maximum 12, as opposed to what happens on the 
torus . In order to allow comparison, we define L eff = (^p) 1//4 , in such a 
way that a hypertorus of L = L e g has the same number of plaquettes as the 
corresponding hypersphere. 
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The specific heat and the Binder cumulant are useful quantities to mon- 
itorize the properties of a phase transition. At the critical point, they are 
known to posses different thermodynamical limits depending on whether the 
transition is first order or higher order, and hence, their behavior with in- 
creasing lattice size can give some clues as to the order of the phase transition. 
We have studied the Finite Size Scaling (FSS) of these two energy cumulants. 

The specific heat is defined for both energies as: 

Cy = = N P ((E 2 P ) - (E P ) 2 ) , (6) 

C2, = ^E 2p = N p ((Ei p ) - <£ 2p > 2 ) . (7) 

As we have observed a high correlation between both energies, the results 
being qualitatively the same, we shall only report on the observable defined 
for the plaquette energy E p . 

In a second order phase transition, scaling theory predicts the specific heat 
maximum to diverge as L a / U . If the transition is first order it is expected 
to diverge like the volume L d (or strictly like the number of plaquettes N p ) 
reflecting that the maximum of the energy fluctuation has the size of the 
volume. This is expected to hold in the asymptotic region of the transition 
L ^> £ c . In the transient region, namely L < £ c , the specific heat is expected 
to grow more slowly than the volume. 

We shall use the number of plaquettes N p to study the specific heat 
behavior in order to have a single parameter for both the HT and the HS 
topologies. This is of course equivalent to doing the discussion as a function 
of L on the HT and L eff on the HS. We have in this case C v ~ (N p ) a / dv . 

We have also studied the behavior of the fourth cumulant of the energy 

VW-g|. (8, 

When the energy distribution describing the system is gaussian Vl — > 
2/3 in the thermodynamic limit. This is the case for a second order phase 
transition. If the transition is first order, far from the critical coupling j3 c , 
Vl also tends to 2/3, reflecting the gaussianity of the energy distribution. 
However, at (3 C the distribution can be described by two gaussians centered 
about the energy of each metastable state E 1 and E 2 , and hence this quantity 
has the non-trivial thermodynamic limit 
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2(Ef + El ] 
3(E? + El 

Another interesting quantity measurable from the density of states is the 



V^oo - 1 - < 2/3 . (9) 



distribution of the partition function zeroes, or Lee- Yang zeroes ||26|| . To 
clarify the role of the partition function zeroes let us first stress a well known 
fact: there are no phase transitions on a finite volume, phase transitions arise 
in the thermodynamical limit and are signaled by non-analyticity in the free 
energy: 

F((3,V) = -^logZ((3) . (10) 

Z(/3) is a linear combination of exponentials, and hence an analytical 
function. This implies that the free energy can be singular only where Z{(3) = 
0. However if the coupling (3 is real, and the volume is finite, that linear 
combination is a sum of positive terms and hence, the zeroes of Z{(3) are 
located in the complex plane of the coupling (3. The onset of the phase 
transition in the limit of infinite volume is signaled by a clustering of the 
zeroes on the real axis at (3 C . 

Here follows the description of our procedure, for a general description of 
the method see [^7j . 

As mentioned previously we work only with E p , and then we consider 
only the spectral density method |28| for this variable. Since we work at 



fixed 7 we consider only the (3 coupling derivatives. 

From the MC simulation at (3 we obtain an approximation to the density 
of states which allows us to compute the normalized energy distribution 
Pp(E). The energy distribution can be expressed as 

Pfs(E) = \w{E)e~? E . (11) 

We use the standard reweighting technique to obtain from the distribution 
measured at (3 the distribution at another coupling lu, which is complex in 
the more general case u = rj + i£. The standard reweighting formula is: 



using ( p]) and the normalization condition Y2e Pp(^) = 1 one obtains 
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fM^P^e-^. (13) 

E 

We can factorize the contributions from the real and imaginary part: 

|M = "£P p (E)(cos(EZ) +i S m(EO)e'^ E . (14) 

This is the standard reweighting formula extended to the complex pa- 
rameter space of couplings. As a first observation we have a pure oscillating 
factor due to Im{uj) 7^ 0. Since E is 0(V) this is a rapidly oscillating function 
which makes it impossible to locate zeroes with large imaginary part. 

The real part of the coupling contributes to the well known exponential 
damping, e~^~^ E ^ which is telling us that we can trust the extrapolation 
only in a small neighborhood (3 ± r\. 

Since Z(Reu) has no zeroes in a finite volume, an easy way to locate the 
zeroes numerically is looking at the minima of the function |G(u;)| 2 , where: 



G{U) ~ W^) ~ E E ME)e-^ • (15) 

The function to minimize is: 



2 (E E P^E)e-^ E cosm) 2 + (E E Pp(E)e-^ E sinm) 2 
' 1 )l (E E P P (E)e-(v-m)2 

(16) 

One interesting property of the partition function zeroes concerns the 
estimation of critical exponents. Denoting by c^o the coupling where the first 
zero is located, the distance to the real axis scales with the v exponent: 

Im(oj ) ~ L- 1/u . (17) 

We shall use this property to compute the effective v exponent. 

From Pe{(3)li we can measure the free energy gap, AF(L), which is the 



difference between the minima and the local maximum of the free energy |24 
We use the spectral density method to get, from the measured histograms, a 
new histogram where both peaks have equal height. We take the logarithm 
of those histograms and measure the energy gap. A growing AF(L) in the 
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Figure 1: Phase diagram of U(l) in the extended (/?, 7) parameter space. The 
dotted lines have been taken from [11]. The symbols correspond to our simulations 
on the torus. The crosses correspond to the peak of C V (L = 16), the errors are 
not visible in this scale. The diamonds have been obtained with hysteresis cycles 
in L = 8. 



asymptotic region of the transition implies a first order phase transition. An 
increase proportional to L d ~ l is expected For the transition to be second 
order, AF(L) must stay constant with increasing lattice sizes. 



3 Numerical Simulation 

Most of the work has been done by simulating the subgroup Z(1024) C U(l), 
since the phase transition associated to the discrete group lies safely far away. 
The overrelax algorithm can be applied only in the simulations with 7 = 
so that the gain in statistical quality due to the overrelax effect is limited to 
the Wilson action, where we have simulated both, the full group U(l) and 
the discrete one for the sake of comparison. For 7 ^ we have simulated the 
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7 = +0.2 


N 


L c s 




r 




P*(L) 


lm(u) Q ) 


4 


5.383 


0.8910 


70 


29142 


0.8906(2) 


0.0169(1) 


5 


7.150 


0.8855 


210 


8000 


0.8855(2) 


0.00604(6) 


6 


8.921 


0.8834 


520 


3800 


0.88330(5) 


0.00244(6) 


7 


10.694 


0.8818 


790 


2300 


0.88182(3) 


0.00090(3) 


8 


12.469 


0.88095 


930 


1930 


0.88097(3) 


0.00040(2) 


7 = 


N 


L c s 




T 


N T 


(3*(L) 


Im(wo) 


6 


8.921 


1.0128 


240 


8300 


1.01340(5) 


0.00616(8) 


8 


12.469 


1.0125 


400 


2300 


1.0127(2) 


0.00196(4) 


10 


16.021 


1.0120 


780 


1200 


1.01212(3) 


0.00075(3) 


12 


19.574 


1.0119 


850 


1150 


1.01194(2) 


0.00031(2) 


14 


23.123 


1.0117 


920 


1900 


1.01168(2) 


0.00013(1) 


7 = -0.2 


N 


L c s 




T 


N T 


(3*(L) 


Im(u;o) 


6 


8.921 


1.1587 


160 


2000 


1.1587(4) 


0.0107(3) 


7 


10.694 


1.1600 


320 


2800 


1.1597(1) 


0.00550(7) 


8 


12.469 


1.1597 


510 


1000 


1.1603(2) 


0.00321(4) 


10 


16.021 


1.1602 


680 


1500 


1.1604(2) 


0.00171(2) 


12 


19.574 


1.1604 


820 


1200 


1.1602(1) 


0.00083(1) 


14 


23.123 


1.1605 


900 


1100 


1.16048(5) 


0.00044(1) 


16 


26.684 


1.1604 


1150 


1200 


1.16038(2) 


0.00023(1) 



Table 1: Statistics of the data obtained for the HS topology 

discrete group since so the simulation is considerably speeded up. 

For every lattice size we consider, we perform trial runs to locate the 
coupling (3*{L) where the specific heat shows a peak. We use the standard 
reweighting techniques to extrapolate in a neighborhood of the simulated 
coupling. Once the peak is located within an error in the fourth digit of 
the coupling, we perform an intensive simulation there to get Pe{/3*)l- The 
statistics performed at these pseudo-critical couplings are reported in Table 
|l] and Table § for the HS and the HT topology respectively. 

Typically we measure the energies every 10 MC sweeps in order to con- 
struct the energy histogram (|i~l~D From this energy distribution we obtain the 
cumulants of the energy we are interested in, and the critical exponents using 
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7 = -0.1 


L 




r 




13* (L) 


Im(uj ) 


6 


1.0720 


350 


1900 


1.0716(2) 


0.0097(1) 


8 


1.0784 


640 


1400 


1.0786(2) 


1.1539(2) 


12 


1.0820 


820 


750 


1.0818(1) 


0.00114(2) 


16 


1.08278 


930 


900 


1.0827(1) 


0.00040(2) 


20 


1.0833 


1150 


1100 


1.0833(1) 


0.00020(1) 


7 = -0.2 


L 


/3sim 


r 


N T 


f3* {L) 


Im(u ) 


6 


1.1460 


380 


2000 


1.1452(2) 


0.0115(1) 


8 


1.1535 


620 


1900 


1.1539(2) 


0.00506(6) 


12 


1.1582 


840 


1200 


1.1582(2) 


0.00152(3) 


16 


1.15935 


920 


900 


1.1593(1) 


0.00060(2) 


20 


1.1599 


1150 


900 


1.1599(1) 


0.00028(1) 


7 = -0.3 


L 


Aim 


r 


N T 


(3*(L) 


Im(u ) 


6 


1.2255 


340 


2000 


1.2237(4) 


0.0138(1) 


8 


1.2344 


560 


1800 


1.2340(1) 


0.00623(4) 


12 


1.2395 


770 


1000 


1.2395(2) 


0.00190(3) 


16 


1.2410 


900 


900 


1.2410(1) 


0.00084(2) 


20 


1.2416 


1100 


1200 


1.24156(5) 


0.00041(2) 


24 


1.2417 


1200 


1100 


1.24162(5) 


0.00022(1) 


7 = -0.4 


L 


/3sim 


r 


N T 


(3*(L) 


Im(u ) 


6 


1.3090 


400 


1800 


1.3082(4) 


0.0152(1) 


8 


1.3192 


600 


1600 


1.3194(3) 


0.00704(5) 


12 


1.3258 


710 


1400 


1.3259(1) 


0.00238(3) 


16 


1.32775 


840 


900 


1.3278(1) 


0.00108(2) 


20 


1.3285 


930 


1600 


1.3284(1) 


0.00054(2) 


24 


1.3286 


1150 


1500 


1.3286(1) 


0.00029(1) 



Table 2: Statistics of the data obtained in the HT topology. 

Finite Size Scaling techniques. We remark that one of the main sources of 
systematic error when measuring critical exponents, is the indetermination 
in the coupling (5*{L) where to measure. 
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The simulations have been done in the RTNN machine consisting of 32 
Pentium Pro 200MHz processors. The total CPU time used is the equivalent 
of 6 Pentium Pro years. 

We have updated using a standard Metropolis algorithm with 2 hits. The 
acceptance has always been between 65% and 75%. 



In order to consider the statistical quality of the simulation, following |29 
we define the unnormalized autocorrelation function for the energy 

1 N-t 

C(t) = 7 =—S2EkE iH -(E)\ (18) 



N-t 

i=i 

as well as the normalized one 

M - W) • (19) 

The integrated autocorrelation time for the energy, r mt , can be measured 
using the window method 



t'=i 



for large enough t, which is in practice selected self-consistently. We use t 
in the range 5r mi , 10r mt , and we check that the obtained r mt remains stable 
as the window in t is increased. 

We have always started from hot and cold configurations in order to make 
sure that the system does not remain in a long living metastable state, which 
could be interpreted as a Dirac sheet. The results coming from both types 
of starts have always been indistinguishable. 



4 The phase transition line $.(7) 

We have studied the deconfinement-confinement phase transition line at sev- 
eral values of 7 (Figure [TJ). 

In the region of negative 7, apart from the deconfinement transition line, 
there is another transition line provoked by the competing interaction be- 
tween the couplings (/3, 7). In the limit 7 = —00 the model is not dynamical 
since for all finite (3, cos# p = 0. So we expect this line to end at the corner 
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Figure 2: MC evolution of E p at 7 = +0.2 on the HS topology for N = 6 (lower 
window), N = 7 (middle) and N = 8 (top). 



(7 = — 00, (3 = +00). We have not performed a deep study of this transi- 
tion line, however the simulations at L = 8 revealed double-peak structures 
pointing to a first order character. The limit 7 = +00 is equivalent to a Z{2) 
theory and the critical point can be calculated exactly by self-duality. 

We have studied carefully the region between lines A and B finding no 
signatures of the existence of additional lines. 

We have focused on the transition line A, at several values of 7, on the 
hyper-torus and on the hyper-sphere. 

The structure of this section is the following: 

First, we describe the results at 7 = +0.2 on the spherical topology. Our 
purpose is to observe how does the HS topology behave in a non controversial 
region, and comparing with the known results on the hypertorus. 

Second, we study the point 7 = (Wilson action) on the spherical lattice. 
This coupling has been recently claimed to be the starting point of a critical 
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Figure 3: E p distribution at 7 = +0.2 on the HS topology. The small window is 
the distribution we obtained on the HT topology in L = 4 at = 0.8595 



line (infinite correlation length at the critical point) which extends in the 
range 7 < 0, with an associated critical exponent v ~ 0.37 fl7l . Since this 
assertion relies on the absence of two-state signals on spherical lattices up to 
N = 10, we shall check whether or not double peak structures set in when 
larger spheres are considered. 

Third, we go to the region of negative 7. We are aware of no previous 
systematic study on the hypertorus in this region, so we run simulations up 
to L = 24 at 7 = -0.1, -0.2, -0.3, -0.4. Motivated by the results on the 
hyper-torus we have just run a single 7 negative value on the hyper-sphere. 
For the sake of comparison with [T7]] we choose this value to be 7 = —0.2. 

Finally, we discuss the Finite Size Scaling behavior exhibited by both 
topologies. 
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Figure 4: Binder cumulant at 7 = +0.2 on the HS topology in N = 4, 5, 6, 7, 8. 
The cross in the upper corner signals the second order value 2/3. 

4.1 Results at 7 = +0.2 on the spherical topology 

As we have pointed out, there is a general agreement on considering the 
phase transition first order in the region of positive 7. 

The lattice sizes used range from N = A to N = 8, which correspond to 
L e g ~ 5 and L c fr ~ 12 respectively see (Table [[]). 

In Figure 0, the MC evolution for N = 6, 7, 8 is shown. We remark that no 
multicanonical update is needed to obtain a very high rate of flip-flops up to 
N = 8. However, on the torus, for L > 6 the probability of tunneling between 
both metastable states is so tiny that a reasonable rate of flip-flops is not 
accessible to ordinary algorithms. Probably the inhomogeneity of the sphere 
is increasing the number of configurations with energies which correspond 
neither to the confined nor to the deconfined phase, but in between. These 
configurations make the free energy gap to decrease and hence the tunneling is 
easier on the HS topology. In short, the sites without maximum connectivity 
act as catalysts of the tunneling. 
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Figure 5: Specific heat maximum, and Ferrenberg-Swendsen extrapolation (solid 
line) at 7 = +0.2 on the HS topology. The small window represents C™ ax (iVp) The 
dotted line corresponds to the behavior expected in a first order phase transition. 



In Figure || the energy distributions are plotted. The distribution at 
N = 6 is distinctly non-gaussian, but a blatant two-peak structure is observed 
only from N = 7 (L e g ~ 11) on. So, when comparing with the result at L = 4 
on the torus, (see small window in Figure |3|) a first observation is that the 
onset of a two-state signal is particularly spoiled by the HS topology, at least 
in 4D pure compact U(l) gauge theory. 

Concerning the latent heat, we remark its stability already at iV = 7, or, 
if anything, its increase from N = 7 to N = 8. One would even say that the 
same happens between the positions of the would-be two states in iV = 6 and 
the position observed in iV = 7. The behavior of the Binder cumulant reflects 
this fact (see Figure |J). We observe a rapid growth of V™ n for small lattices 
sizes, apparently towards 2/3 (second order value). This growth stops when 
the splitting of the two peaks is observed. The splitting of the two peaks is 
reflected by smaller values in Vf am . We remark that the errors quoted for 
ymm are ca i cu i a ted taking into account the indetermination in the value of 
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the minimum, but not the displacement in the position of the coupling where 
the minimum appears. 

In our opinion, mainly two reasons can give account of this behavior. 
The first one comes from general grounds: at the very asymptotic region 
of a first order phase transition, the energy jump gets sharper and sharper, 
and a slightly increase of the latent heat could be expected. The second one 
would be the increasing restoration of homogeneity in the hypersphere when 
increasing the lattice size. The latent heat observed in small N might be 
affected by the inhomogeneity of the hyper-sphere. 

At this point we cannot give a single reason for this to happen. We 
postpone a stronger conclusion to the section devoted to the Finite Size 
Scaling discussion. 

The value of the latent heat when obtained with a cubic spline fit to the 
peaks in N = 8 is Ci at = 0.064(2). Results obtained with mixed hot-cold 
starts in N = 9 seems to give an energy jump around 0.067. However the 
reliability of this method is very limited and we do not dare to extract strong 
conclusions from it. Taking into account that he cubic spline at N = 7 gives 
0.053(2) a possible scenario would be a slowly increasing latent heat towards 
its asymptotic value Ci at (oo). 

The peak of the specific heat for the different lattice sizes is displayed in 
Figure [5[ The continuous line represents the FS extrapolation. The plot of 
the peak value, C^ 1SX (N P ) as a function of the plaquette number reveals a 
linear relationship, and hence a first order character. 

4.2 Results at 7 = on the spherical topology 

Let us first situate the status of the studies using the toroidal topology with 
the Wilson action. As we have pointed out in the introduction, the tran- 
sition was believed to be continuous till the simulation of a L = 6 lattice 
revealed the first signs of the existence of two metastable states H. Numeri- 
cal simulations up to L = 10 showed a rapid decrease of the latent heat when 
increasing the lattice size, suggesting that the two metastable states might 
superimpose in the limit L — > 00 ||. However, in |30|| , Azcoiti et al. suggest 
that the latent heat starts to stabilize at L ~ 12, 14. 

When one tries to simulate large lattices the tunneling becomes scarce 
and the use of multicanonical simulations for those lattice sizes seems to be 
in order. However, in spite of its usefulness in spin models, we are aware of no 
multicanonical simulation showing a substancial flip-flop rate improvement 
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for this model. An alternative procedure to enhance the tunneling probability 
has been proposed in p2| . Within this approach the use of a monopole term 
in the action as a dynamical variable, has been revealed as a useful method 
to increase the flip-flop rate. 

Lattice sizes up to L = 16 have also been studied using multihistogram- 
ming techniques || and RG approaches || [Kj . The results support the idea 
of a quasi stable latent heat for L > 12. 

Topological considerations stressed already in the introduction led some 
authors to use lattices homotopic to the sphere. In the case of the Wilson 
action (7 = 0) the two-state signal is absent up to iV = 10 [fT7| . Following 
this, we have simulated on the spherical topology in lattices ranging from 
N = 6 to iV = 14, finding that the two-state signal sets in from N = 12 on. 
The preliminary results up to iV = 12 are already published f23 |. We include 
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Figure 6: MC evolution of E p at 7 = on the HS topology for N = 10 (lower 
window), N = 12 (middle) and N = 14 (top). The vertical dashed lines separate 
different runs 
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Figure 7: E p distribution at 7 = on the HS topology. The small window 
corresponds to our simulation in the HT topology in L=6 at (3 = 1.0020. 



here the results with higher statistics in iV = 12 and the results in N = 14. 

The MC evolution for iV = 10, 12, 14 is plotted in Figure |[ For the 
N = 14 lattice we have four independent runs, signaled by the dashed lines 
in the figure, all of them giving the same predictions. 

In Figure [7| the distribution of E p is plotted in lattices ranging from N = 6 
to N = 14. We observe that a two-peak structure is revealed first time by the 
histogram in N = 12, which has a L c g ~ 19. On the toroidal topology the 
equivalent signal is observed already at L = 6 (see small window in Figure 
0). 

Having in mind the results of the previous section at 7 = +0.2, it seems 
that the minimum L e g- required to observe a two-peak structure in the spher- 
ical topology is around three times the minimum lattice size L needed in the 
torus to observe two peaks. 

The lattice sizes used in |l7j at 7 = ranged from N = 4 to iV = 10, so 
it is not surprising that they did not see any two-peak structure. Our results 
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Figure 8: E p distribution at 7 = on the HS topology simulating the full group 
U(l). 



are in that sense compatible with theirs, though ours show a faster divergence 
for C™ ax (iV) with increasing N than the one observed in [17] in the lattice 
sizes we share (N=6,8,10). However we are particularly confident on this 
point because our simulation has been performed closer to the peak of the 
specific heat P*(L), and so we expect the Ferrenberg-Swendsen extrapolation 
to be more precise. 

For the sake of comparison with the full group U(l), simulations with the 
Wilson action have been performed at 7 = 0. At this coupling an imple- 
mentation of overrelax is possible, and the global decorrelating effect should 
manifest itself in a better statistical quality. 

In Figure |8| we show the histograms of E v from iV = 6toiV = 12 sim- 
ulating the full group. The results are fully compatible with those obtained 
from simulations with Z(1024). 

This test being performed, we go back to the description of the results 
obtained for the discrete group. 
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Figure 9: Specific heat maximum, and Ferrenberg-Swendsen extrapolation (solid 
line) at 7 = on the HS topology. The small window represents C™ ax (iVp) The 
dotted line corresponds to the behavior expected in a first order phase transition. 



In Figure |9| we plot C V (N) for different N's, together with the Ferrenberg- 
Swendsen extrapolation. The small window is a log — log plot of C™ ax (A/') 
as a function of N p . A linear behavior is observed from N = 10 on, which 
is even faster than linear when iV = 14 is taken into account. This fact by 
itself implies the first order character of the transition since it means that the 
maximum in the energy fluctuation has the size of the volume. As a further 
check we have measured the v exponent from the scaling of the Fisher zeroes 
(see Table ||). 

In Figure |10| we plot \G(u) | 2 for several lattice sizes on the spherical topol- 
ogy at 7 = 0. The different curves stand for the different uj we extrapolate. 
For small lattice sizes Im(u)o) is larger, and the damping is more severe than 
for the larger lattices since in the later the imaginary part contributes with 
a faster oscillating function. Actually, we observe that for N = 14 even a 
second minimum could be measured before the signal is damped, while for 
N = 10 one can measure accurately only the first one. 



22 



^3 



0.5 z 

0.4 E- 

0.3 — 

0.3 — 

0.1 — 

0.0 = — 
0.0000 
0.3 



0.2 
0.1 



0.0 
0.0000 
0.020 



0.015 
0.010 
0.005 
0.000 



0.0006 



.„ii... 



0.0002 



IliiiUy 

0.0004 



0.0002 



0.0004 



0.0006 



Mill 



"HUlll 



0.0008 



0.0010 



Im(c 



illini um 

0.0006 





1 1 


1 1 1 1 1 1 


Ji 


ilil 

. ill' 11 " 


i " i1 ! iiii1i 'iiiiiinii II 1 - 



0.0008 



1 14 



0.0012 



Figure 10: Plot of (jl|) for N = 10 (lower window), N = 12 (middle) and N = 14 
(top) at 7 = in the HS topology. 



The results for the effective v are plotted in Figure |TT] lower window. As 
could be expected from the behavior of the specific heat, for small lattices the 
effective exponent is somewhat larger than 0.25. It gets compatible with the 
first order value from A" = 10 on. In the largest lattices the distance between 
the two peaks slightly increases, and we measure a u e s slightly smaller than 
the first order value. As a v < 0.25 is strictly impossible we expect this to be 
a transient effect due to finite size effects associated to the observed splitting 
up of the two peaks. 

Concerning the latent heat, we observe a behavior completely analogous 
to the one observed at 7 = +0.2. We do observe a two peak structure quite 
stable when comparing AT = 12 and A" = 14, but it slightly increases in 
= 14. The plot of the Binder cumulant reflects again this fact (see Figure 
|i~2|) . The fast growth of V™ n towards 2/3 is preempted by the onset of double 
peaked histograms from = 12 on, and it even decreases in A^ = 14. 
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Figure 11: Effective exponent v at 7 = (lower window), and at 7 = —0.2 (upper 
window) on the spherical topology. In the smaller lattices an v e g ~ 1/3 is observed 
which becomes 1/d when large enough lattices are considered. 



The cubic spline fit in N = 14 at the peaks gives for the latent heat C\ & t ~ 
0.018(2). which is compatible with the results suggested by extrapolating to 
infinite volume the values obtained on the torus up to L = 14. Our results are 
hence supporting the conjecture stressed in |3(J about a quasi stabilization 
of the latent heat on the torus from L — 12 on. 



4.3 Results for 7 < 
4.3.1 Toroidal topology 

We have performed a systematic study of the transition at several negative 
values of 7 on the toroidal topology. 
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Figure 12: Binder cumulant at 7 = (Wilson action) on the HS topology in 
N = 7, 8, 10, 12, 14. The cross in the upper corner signals the second order value 



2/3. 



We find that the two-state signal persists for all 7 values we consider. 
The energy histograms reveal an increasing weakness of the transition when 



going to more negative 7 values (see Figure [L3|). A double peak structure 
is clearly visible at 7 = —0.1 in L = 8, while at 7 = —0.4 one has to go to 
L = 12 to observe an equivalent signal. 

In what concerns the behavior of the free energy gap AF(L), it grows 
for all investigated lattice sizes at all negative 7 values (see Figure [14] lower 
window). The value of L at which AF(L) starts growing is certainly larger 
as the value of 7 is more negative. This is another test of the increasing 
weakness of the transition as 7 gets more negative. 

The statistics performed on the torus are reported in Table [| We also 
quote for the different negative 7, the value of (3*{L) and the position of 
partition function zero closest to the real axis. We have computed from the 
imaginary part of the zeroes the effective v exponent between consecutive 
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Figure 13: E p distributions at the different negative 7 on the HT topology. We 
have run up to L = 20 at 7 = -0.1, -0.2 and up to L = 24 at 7 = -0.3, -0.4. 



lattice sizes following ([T7|). 

In Figure |T3| we plot for the different 7 values the we measure. In 
all cases a v e ^ ~ 1/3 is observed for small lattice sizes, which gets closer to 
0.25 when the lattice is large enough. From this figure the trend of v e s seems 
rather clear towards the first order value. 

From the energy distributions we measure the latent heat through a cubic 
spline at the peaks. Taking into account the value z/ e fj = 0.25 we measure, 
we plotted it as a function of the inverse of the volume L~ 4 , which is also the 
expected behavior of the latent heat when the transition is first order. The 
latent heat can be extrapolated to a value which is safely far from zero (see 
Figure |TJ| upper window). 

In Figure [16] we plot C™ ax (L) as a function of N p . As could be expected 
from the behavior of the effective exponent z/, the maximum of the specific 
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Figure 14: Free energy gap AF(L) (lower window) and latent heat (upper win- 
dow) for the different negative 7 in the HT topology. 



heat for small lattices diverges slower than the volume. For the smaller lat- 
tices the effective exponent is ajv ~ 1.4. This value increases monotonically 
with the lattice size. In the largest ones we observe otjv ~ 3.5. 

4.3.2 Spherical topology 

On the toroidal topology we have found that the minimum lattice size re- 
quired to observe a two state signal is obtained through an apropriate com- 
bination (7, L m in), with increasing L min for increasingly negative 7, the be- 
havior being qualitatively similar for all the 7 values we have investigated. 
In view of this, we have studied a single 7 < value on the spherical lattice 
to check if the two state signal is absent on this topology. For the sake of 
comparison with WB, |I7[ we choose this value to be 7 = —0.2. 



We have run simulations on spheres ranging from N = Q (L e fj ~ 8) to 
N = 1Q (L eS ~ 26). The MC evolution is plotted in Figure [T7|. 

The distribution in iV = 14 is distinctly non gaussian and the splitting 



of the peaks occurs in iV = 16 (see Figure 18). We remark that the sim- 



ulations in N = 16 are extremely expensive in CPU. In order to alleviate 
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Figure 15: Effective ^ exponent for the different negative 7 in the AT topology. 



thermalization we have parallelized the code using shared memory in two 
PPro processors. In this lattice we have run two independent simulations 
starting from different configurations the results being fully compatible. 

The behavior of the maximum of the specific heat is shown in Figure [19|. 
We observe the same trend than in the previous values of 7 in the sphere, 
i.e., an increasingly fast divergence of the specific heat with the lattice size, 
showing an effective ajv ~ 4 already when the histogram width becomes 
constant with increasing lattice size. 

At this value of 7 we have a worse estimation for the latent heat to be, 
as the two peaks have not split enough to allow an accurate measurement. 

Again the behavior of the Binder cumulant is very significant (see Figure 
2D[ The value V™ n shows a very fast trend towards 2/3 in the small lattices. 
When increasing the lattice size the rate gets slower, and finally the value in 
N = 16 is compatible with the value in N = 14 preempting the extrapolation 
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Figure 16: C™ ax (L) as a function of -/V p for the different negative 7 values on the 
HT topology. We have used log — log scale for the clarity of the graphic sake. The 
dotted line corresponds to the first order behavior. 



to 2/3. Unfortunately lattices larger than N = 16 are unaccessible to our 
computers nowadays, and we cannot observe a decreasing V™ n as we did for 
the other 7 values. However from the behavior between iV = 14 and iV = 16 
an increasing V™ m for larger spheres seems to be very unlikely. 

Concerning the effective critical exponent v we have measured the posi- 
tion of the first Fisher zero (see Table [l|) and compute (see Figure pT|) . 
Again the first order value 1/d is reached in the largest lattices. 



4.4 Toroidal versus Spherical topology 

Due to the fact that on the sphere there are a number of points with less 
than the maximum connectivity, violations to standard FSS in the form of 
uncontrolled finite size effects are expected to appear on this topology. 

We have observed all along this work that working on the spherical topol- 
ogy retards the onset of two-state signals. 
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Figure 17: MC evolution of E p in N = 12 (lower window), N = 14 (middle 
window) and N = 16 (upper window) at 7 = —0.2 on the HS topology. The two 
different runs in N = 16 are separated by a vertical dashed line. 



In terms of L e g, some qualitative prediction arising from our results would 
be that at fixed 7, the minimum lattice size required to observe a double peak 
structure, L min , is around three times larger on the sphere than on the torus. 
This one being the showier difference, is not certainly the only one. 

When comparing C™ ax (L e g) on both topologies one finds always smaller 
values on the spherical topology (see small window in Figure ^l]). We did 
run on L = 12 on the torus and iV = 8 on the sphere, which has an L c g ~ 12, 
sweeping an interval of (3 values including the phase transition. In Figure ^1 
we have plotted E p in both topologies. In the region of low (3 the system is 
disordered, the entropy is higher and the system is not so sensitive to inho- 
mogeneities of the lattice. However, in the region of high f3 the energy is in 
general smaller than for the homogeneous lattice since the system tends to be 
ordered, and the influence of the sites with less than maximum connectivity 
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Figure 18: E p distributions at 7 = -0.2 on the HS topology. 



is more evident. This difference gets smaller when the lattice size is increased 
on the sphere, and hence the energy jump is larger. This feature explains the 
observed splitting of the two states on the spherical lattices when increasing 
lattice size. 

Altogether, it is not recommended to work on spherical lattices to check 
the existence or not of two-state signals. However, there are a number of facts 
that make this topology not so dissapointing in spite of those uncontrolled 
finite size effects. The first one concerns the behavior of (3*(L) and the second 
one the measure of the latent heat. 

In Figure ^2] we compare (3*{L) on the torus and on the sphere. A first 
observation is that despite they not having the same values in finite lattices, 
both curves get closer when increasing the lattice size. This supports the 
idea of a common thermodynamic limit for both topologies. On the other 
hand, the shift in the apparent critical coupling with L is much less dramatic 
on the sphere than on the torus. It seems that finite size corrections to f3 c 
are smaller in the case of the spherical topology. This behavior has also been 
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Figure 19: Specific heat maximum, and Ferrenberg-Swendsen extrapolation (solid 
line) at 7 = —0.2. The small window represents C™ ax (7V p ) The dotted line corre- 
spond to the behavior expected in a first order phase transition. 



observed in the Ising model on spherical lattices [31 



The behavior of the latent heat is very significant. In Figure 23 



we com- 
pare the latent heat in finite volumes for both topologies. On the Wilson line, 
at 7 = the asymptotic value of the latent heat is obtained in the spherical 
topology. The inhomogeneity of the sphere has, paradoxically, helped us to 
run lattices with L e s ~ 24 without having to worry about the tunneling rate. 
From this graph we quote Ci at (oo) = 0.018(1). In the upper window the 
same is plotted for 7 = —0.2. Unfortunately we have just a single lattice 
size on the sphere to measure the latent heat, however, from the behavior 
exhibited by this topology, an increase of the latent heat on spheres larger 
than N = 16 seems rather likely. Altogether, at 7 = —0.2 a Ci at (oo) ~ 0.009 
is plausible. 

In what concerns the latent heat, the spherical topology seems to afford 
an useful complement to the results obtained for the the toroidal lattice. In 
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fact, on the Wilson line lattices larger than L — 16 cannot be studied due 
to the technical problem associated to the scarce tunneling. The spherical 
lattice alleviates this technical problem, and the results are supporting the 
value of the extrapolated latent heat as a function of the inverse of the 
volume, from the data on toroidal lattices up to L = 16. One could consider 
the possibility of simulating on spherical lattices at negative 7 to solve the 
tunneling problem, but in our opinion, the price to pay is too high. 

5 Conclusions 

The first order character of the deconfinement transition in pure U(l) has 
been proved, up to the limits of a rather reasonable numerical evidence, in 
the interval 7 G [+0.2, -0.4]. 

In 7 = +0.2 we have been able to stabilize the latent heat. We are aware 
of no simulation on the torus showing a stable latent heat due to the scarce 
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Figure 21: £" p at 7 = -0.2 in L cff ~ 12. The dotted line corresponds to the 
toroidal topology, the dashed one to the sphere N = 8. The small window shows 
the difference in the specific heat between both cases. 



tunneling in lattices larger than L — 6. The spherical topology has helped 
to solve this problem. However, probably any lattice with inhomogeneities 
would produce the same catalysing effects. 

In 7 = we have also been able to measure the asymptotic value of 
the latent heat. We have proved the suggestions of several authors about a 
quasi stabilization on lattices larger than L ~ 12. The data on the spherical 
topology have been crucial to discard the possibility of a slowly vanishing 
latent heat. It follows that the discretization of pure compact U(l) LGT on 
the lattice using the Wilson action exhibits a first order phase transition with 
a latent heat in the thermodynamical limit Cut ~ 0.018. 

On the toroidal topology things happen qualitatively in the same way 
than in 7 = up to 7 = —0.4. We have run on spherical lattices in 7 = —0.2 
looking for an argument similar to the one found in the Wilson case concern- 
ing the stabilization of the latent heat. The simulation in N = 16 gives an 
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Figure 22: f3*(L) for the different lattice sizes on the torus (dotted lines) and 
on the sphere (dashed lines) at 7 = (lower window) and at 7 = —0.2 (upper 
window). The couplings for 7 = on the torus have been taken from fflOfj . 



estimation for the latent heat though rather imprecise because the splitting 
of the peaks is not good enough for the measurement to be accurate. How- 
ever, in view of the behavior exhibited by the spherical topology concerning 
the trend of the two peaks on spherical lattices to separate, we are prone to 
consider the value measured in iV = 16 as a lower bound for the latent heat 
in 7 = —0.2. On the other hand, we would at present say that the latent heat 
extrapolated from the data on the torus up to L = 20 is rather accurate, in 
view of the behavior for the Wilson case. The possibility of running larger 
spheres surpasses our computer resources. A highly parallelized version of 
the code should be used to alleviate thermalization, which is a possibility we 
do not discard completely at medium term. 
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Figure 23: Latent heat for the different lattice sizes on the torus and on the 
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Either proving or discarding the possibility of a TCP at some finite neg- 
ative |7tcp| > 0.4 will be a very difficult task from the numerical point of 
view. An analytical argument would be welcome. A very tiny two-state sig- 
nal, comparable with the one observed in L = 8 at 7 = —0.4, is observed 
on the torus in L = 16 at 7 = —0.8, which is the most negative value we 
have run on the torus. Taking into account the factor 3 in L e ^ one would 
need a iV ~ 30 sphere to observe a tiny double peak structure. Simulation 
with spherical lattices in this range of 7 have no sense, and nothing can be 
concluded from the absence of two state signals from such negative 7 values. 
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In view of these difficulties to stabilize the latent heat for very negative 
7 values, the only chance to discern the order of the phase transition is the 
study of effective critical exponents. Cumulants of the energy, such as the 
Binder cumulant, have been shown to behave like expected in the first order 
case only when the stable two-peak structure is almost setting in, and there, 
we do not need further evidences any longer. We could not expect additional 
information since by definition the Binder cumulant relies on the existence 
of stable latent heat in order to extrapolate to < 2/3. 

The advantage of studying the effective exponents (which is nothing but 
studying the evolution of the histograms width) is that we do not need a 
direct observation of latent heat to conclude that a transition is first order 
3"E|, |3"B|] . We have observed a v e ff which evolves monotonically until it reaches 



the first order value, namely 0.25, in all cases. The statistics needed to 
observe a monotonous behavior, are order lOOOr at the coupling /3*(L), which 
has to be located with high precision (four digits in our experience) in order 
to accurately measure effective exponents. However, one has to make every 
effort to observe the trend of the effective exponents, since it is the only 
chance to discern the order of such tricky transitions. 

In pure U(l), for the Wilson case, an exponent v ~ 1/3 was widely 
predicted 0, |3|, |ll|] at the beginning of the eighties when the lattice sizes 
where too small to reveal two peaks. That v was shown to become v ~ 0.29 
when simulating L = 14 p0| |. We have measured v = 0.25 on spherical 



lattices and stated its first order character. Within the approximation of 
effective potentials it can be shown that along the transient region of a weak 
first order phase transition, everything goes like in a second order one with a 
thermal exponent a = 0.5 f2"0|| . Together with Josephson law (a = 2 — ud) it 
implies v ~ 0.37. This statement has been checked in 2D Potts f20fl , 3D and 
4D O(N) models pi B3, M, 133 and in the 4D SU(2)-Higgs at T=0 with 



fixed Higgs modulus f36fl . Our results prove that pure compact U(l) theory 
behaves in the same way. 

A general remark to draw from this paper is though there is no theo- 
rem constraining transitions presenting z/ e g ~ 1/3 to be of WFO, yet the 
occurrence of such exponent is a warning to consider the phase transition as 
suspicious of being first order. 
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Appendix: Schwinger-Dyson Equations 



As a further check we have implemented the Schwinger-Dyson equations 
(SDE) [^8| which allow one to recover the simulated couplings from the Mon- 
tecarlo data. 

Let A(6) be an operator with null expectation value: 

(A(9)) = Z- 1 [ [d6]A(8)e- s[e] = . (21) 



Derivating with respect to 9 this equation trivially yields 

<T> = (22 > 

which is the equation of movement of the operator A(0) or Schwinger- 
Dyson equation. When the action depends on several couplings the equation 
([22] ) can be expressed as: 

.dA{6) 



i 

This equation relates the values of the couplings with the expectation 
values we measure from the MC simulation. We need as many independent 
equations as couplings we have to determine in the action. In our case, in 
order to measure both j3 and 7, we need two operators with null expectation 
value in order to have two independent tests. At each lattice site n and for 
every direction /1 we consider the operators: 

A(9) = sm6 p = sin(# niM - 6 stap ) , (24) 
B(6) = sin 26 p = sm2(6 rhP - 6 stap ) , (25) 

where 9 sta p is the staple of the link labeled by (n,/x). 
Applying equation fl23|) to those operators we get: 



cos(6> njA1 - Ostap)) = fiC/, sm(9 n ^ - 6 sta p) sin(6*„ iM - 9 sta p)} + 
p p p 

+2 7 (J2 sm(0„, M - 9 stap ) 2(0^ - 6 stap )) (26) 



41 



2 COS 2 (°n,fi ~ Ogtap)) = Sill2 ( d n,fi ~ ^stap) ^ sin(0 n>jU - 6 stap )) + 

P P P 

+2 7 (^ sin 2{9 n ^ - 9 stap ) £ sin 2{6 n ^ - 6 stap )) (27) 
p p 

On the hypertorus ^ p means the sum over the plaquettes in positive and 
negative directions (±/i) bordering the link (12 plaquettes). On the sphere one 
has to be careful, since not all the links have 12 surrounding plaquettes, and the 
sum has to be understood as extended to the existing plaquettes. 

These equations hold for all n in such a way that we can sum up the equations 
for every single link and average over the number of links. We quote in Table [| for 
the tests done at 7 = —0.2 using both topologies. They show a perfect agreement 
between the simulated couplings and the ones recovered from the MC expectation 
values. 



HT topology 


HS topology 


L 




(At)sd 


N 


/3sim 


(/?,7)sd 


6 


1.1460 


1.1466(44),-0.2009(26) 


6 


1.1587 


1.1592(14),-0.1998(12) 


8 


1.1535 


1.1527(22) -0.1997(15) 


8 


1.1600 


1.1599(12) -0.1999(8) 


12 


1.1582 


1.1581(12),-0.1999(8) 


10 


1.1602 


1.1602(7) -0.2001(9) 



Table 3: Couplings obtained from the MC simulations at 7=-0.2 using the 
Schwinger-Dyson equations. 
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